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If two ore more bodies are immersed in a critical fluid critical fluctuations 
of the order parameter generate long ranged forces between these bodies. Due 
to the underlying mechanism these forces are close analogues of the well known 
Casimir forces in electromagnetism. For the special case of a binary liquid 
mixture near its critical demixing transition confined to a simple parallel plate 
geometry it is shown that the corresponding critical Casimir forces can be of 
the same order of magnitude as the dispersion (van der Waals) forces between 
the plates. In wetting experiments or by direct measurements with an atomic 
force microscope the resulting modification of the usual dispersion forces in 
the critical regime should therefore be easily detectable. Analytical estimates 
for the Casimir amplitudes A in d = 4 — e are compared with corresponding 
Monte-Carlo results in d = 3 and their quantitative effect on the thickness of 
critical wetting layers and on force measurements is discussed. 
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I. INTRODUCTION 

The phase diagram of a fluid is influenced by the presence of a surface in many different 
ways. Most prominent is the modification of the critical behavior of a fluid near a wall 
and the occurrence of new phase transitions induced by the wall such as wetting and drying 
0. For binary liquid mixtures external walls usually manifest themselves by a preferential 
affinity of the wall material for one of the components M which in the vicinity of the critical 
demixing point leads to the phenomenon of critical adsorption of the preferred component 
|5|||]. If the system is made finite by the introduction of a second wall or by confining the 
system to another finite geometry the critical behavior of the fluid is modified again if the 
correlation length becomes comparable to the system size where the size dependence 

of thermodynamic functions takes a scaling form. A finite geometry may also by generated 
spontaneously by a critical fluid if, e.g., a binary liquid mixture near its critical demixing 
transition forms a macroscopic wetting layer on the surface of a substrate [0,0. With the 
introduction of the second surface the variety of phenomena in the confined fluid goes far 
beyond critical finite-size scaling. Apart from the shift of the critical point of the system 



|Tll , |12l one encounters the phenomenon of capillary condensation [|13j if the confining walls 
of the film geometry consist of the same material. The confinement of the fluid causes 
the liquid vapor coexistence line to be shifted away from the coexistence line of the bulk 



fluid into the one-phase regime of, e.g., the bulk vapor |12|,[13| . For not too small wall 
separations a first order phase transition occurs from a confined vapor to a confined fluid as 
the undersaturation of the vapor is lowered at fixed temperature. In a constant-temperature 
plane of the phase diagram the line of two-phase coexistence is terminated by a capillary 
critical point characterized by a critical undersaturation and a critical wall separation beyond 
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which capillary condensation no longer occurs |T3[. Fluid layers growing on the inner walls 
of the capillary reduce its effective width and therefore generate correction terms to the 
well known Kelvin equation, which describes the aforementioned shift of the liquid vapor 
coexistence line as a function of the width of the capillary fl4| . [T5 



From the theoretical point of view these phenomena can be described using density func- 
tional theory |13] and computer simulations of lattice gas models [0,0 . These lattice gases 
are equivalent to Ising models, where the presence of the walls is described by surface fields 
which impose a finite surface magnetization on the Ising system. Density or concentration 
profiles of confined fluids or binary liquid mixtures, respectively, then translate to the mag- 
netization profile of the Ising model. For the description of capillary condensation an Ising 
model with surface fields of the same sign is appropiate. The behavior of the system changes 
drastically, if opposing surface fields are considered. For a confined binary liquid mixture 
this means that the walls perfer different components. It turns out that in this case new 
quasi wetting transitions occur which can be first-order, critical, and tricritical and converge 
to the usual wetting transitions for growing wall separation [16]. Furthermore, two phase 



coexistence becomes restricted to temperatures located below the wetting temperature, if the 
surface fields are equal in opposite [|17| , |18| . The scaling behavior of the magnetization profile 



of an Ising model with opposing surface fields and the dependence of the interface position 



on the strength of the surface fields and the temperature has been studied thoroughly [19,20 



Capillary condenstion does no longer occur, instead one observes the interface derealization 
transition, i.e., the interface in the magnetization profile detaches from one of the walls and 
moves to the midplane of the film. This transition is second order and its critical point 
can be identified with the shifted critical point of the confined system, which in this case 
is located on the temperature axis [18]. Above the critical temperature the magnetization 



profiles become perfectly antisymmetric about the midplane of the film. By increasing the 
strength of the surface fields the critical temperature diminishes and in the limit of infinitely 
strong opposing surface fields the interface derealization transition becomes suppressed at 
all. 



Confined critical fluids also generate long-ranged forces between the confining walls [21 



a phenomenon, which is a direct analogue of the well-known Casimir effect in electromag- 
Contrary to the usual dispersion forces, which are still under investigation for 

|2~6|J27|], critical 



netism 22 



bodies with curved surfaces p3^p5| and in presence of surface roughness 
Casimir forces are governed by universal scaling functions P,|28[|. At the bulk critical point 
these scaling functions reduce to the universal Casimir amplitudes [0,^SJ- Especially for the 
strip geometry a variety of exact results are known from conformal invariance p9| . Away 
from the critical point the scaling functions are only known exactly for an Ising model con- 
fined to a strip in d — 2 [Q . In higher dimensions only the spherical model has given access 
to further exact results for the scaling function of the Casimir force ||3l| . For the O(N) 
universality class in d = 3 so far only approximate results are known based on real space 
renormalitation [j32], the field theoretic renormalization group [p8], and Monte-Carlo simu- 
lations for the film geometry. More recently Casimir forces between spherical particles 
immersed in a critical O(N) symmetric systems have been investigated by field-theoretic 
methods augmented by conformal invariance considerations ||34j| . The field-theoretic treat- 
ment of critical systems confined to finite geometries is notoriously difficult, because the 
theory has to interpolate properly between critical behavior in different dimensions. There 
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has been remarkable progress in devising alternative renormalization prescriptions beyond 



the standard minimal subtraction scheme [35| and in contracting effective actions for the 



Ising |36| and the more general O(N) universality class |3~T |. However, these approaches 
have been devised for finite systems with symmetry conserving boundary conditions, their 
implementation for systems with symmetry breaking boundary conditions (surface fields), 
in which we are interested here, is still lacking. Within the framework of Ginzburg-Landau 
descriptions of critical finite systems in presence of surface fields the theoretical treatment 



has been limited to mean-field considerations for the film geometry |g,|IT|,^,^,^] and con- 
centric spheres [[HJ which can be mapped onto two-sphere and wall-sphere geometries at 
the bulk critical point by conformal transformations (see Ref. ||34|). For tri critical systems 



between parallel plates a thorough mean- field analysis has also been performed . In this 



paper we will concentrate on the Casimir forces in critical films in presence of surface fields, 
which is the adequate description for confined binary liquid mixtures |4[]. 

The remainder of the presentation is planned as follows. In Sec. II we introduce the 
field-theoretic model of a confined binary liquid mixture close to its critical demixing point 
and an adequate Ising model for which the Monte-Carlo simulations of the Casimir force are 
performed. Sec. Ill is devoted to a survey of mean-field results for the scaling functions of 
the Casimir force. In Sec. IV we present one-loop results and Monte-Carlo estimates for the 
universal Casimir amplitudes which characterize the strength of the Casimir forces at bulk 
criticality. We restrict ourselves to the bulk critical point, because the one-loop calculations 
are based on the standard e-expansion which cannot cope with the dimensional crossover. 
In Sec. V we discuss implications of the results presented in Sec. IV for force measurements 
and wetting experiments with critical binary liquid mixtures and we summarize the main 
results in Sec. VI. The one-loop calculation requires the knowledge of the mean-field order 
parameter profiles which are rederived and discussed in Appendix A. The eigenmode spectra 
are derived in Appendix B and the regularization of the one-loop mode sums is described in 
Appendix C. 



II. MODEL 

For the analytical part of the current investigation the standard 4 Ginzburg-Landau 
Hamiltonian 7i = 7ifr + H s for a O(N) symmetric critical system in a parallel plate geometry 
is used. Specifically, the model is defined by the bulk Hamiltonian 

n b m = jd d - i xj Q L dz[ 1 -(v^r+^+^) 2 }, (2.i) 

where L is the Film thickness, <fr = (3>i(x, z), . . . , $at(x, z)) is the N component order 
parameter at the lateral position x and the perpendicular position z (0 < z < L), r is the 
bare reduced temperature, and g is the bare coupling constant. The presence of the surfaces 
gives rise to the surface contribution 

= J d d ~ l x {f [*(x, 0)] 2 + |[$(x, L)f - h x ■ $(x, 0) - h 2 ■ <&(x, L), } . (2.2) 

to the Ginzburg-Landau Hamiltonian, where c\ and c 2 are the surface enhancements which 
characterize the surface universality class J|. In mean field theory and within the dimen- 
sional regularization scheme for the field-theoretic renormalization group q > defines the 



3 



ordinary (O) surface universality class and q < defines the extraordinary (E) surface uni- 
versality class. The leading critical behavior of a semiinfinite system with an O or an E 
surface is described by the two stable renormalization group fixed point values c = +00 and 
c = —00, respectively. Finite positive or negative values of Cj only yield corrections to the 
leading behavior. Within this setting c = is an unstable fixed point, so that (r, c) = (0, 0) 
has meaning of a multicritical point at which both the bulk and the surface of a semiinfi- 
nite system simultaneously undergo a second order phase transition M. This mulitcritical 
point defines a surface universality class in its own right which is commonly denoted as the 
surface-bulk (SB) or special universality class. In the language of a spin model c denotes 
the deviation of the exchange interaction between spins in the surface from its value at the 
multicritical point [see also Eq. (|2.3| ) below]. 

The quantities hi and Ii2 denote surface fields which explicitly break the O(N) symmetry 
of the model. In case of a broken symmetry at the surface in principle also cubic surface 
fields need to be considered [[J . However, for the investigation of the leading critical behavior 
in the presence of nonzero linear surface fields cubic surface fields can be disregarded |J . 

As pointed out in Sec. I a wall which is in contact with a binary liquid mixture will in 
general show some preferential affinity for one of the components so that the composition 
profile varies as a function of the perpendicular coordinate z. This situation can be repre- 
sented by setting c\ > and c 2 > in Eq. ( j2.2|) and prescribing finite values for the surface 
fields hi and h 2 . The phase transition in the bulk in presence of nonzero surface fields is 
called the normal transition f42|| . As far as the leading critical behavior is concerned the 



normal transition is equivalent to the usual extraordinary transition PJ42"|, which can be 



represented by setting hi = h 2 = and choosing c\ < and c 2 < 0. In the following we will 
therefore exclusively use the surface field picture of the extraordinary transition. 

In the field-theoretic analysis only the cases of strictly parallel and strictly antiparallel 
surface fields hj = (hi, 0, . . . , 0), i — 1, 2 will be considered. For the leading critical behavior 
it is sufficient to discuss only the limiting cases hi,h 2 — > ±00 |2j. The above restriction 
to parallel and antiparallel surface fields then means that we only consider the two cases 
hi — hi — > +00 and hi = —h 2 — > +00. To simplify the notation we will refer to the 
former the (+, +) boundary condition and to the latter CclSG clS the (+, — ) boundary 

condition which are the only combinations (E, E) of the E surface universality class in the 
film geometry considered here. One can also combine a symmetry breaking E surface with 
a symmetry conserving O or SB surface. However, as will be demonstrated below, the 
combinations (O, E) and (SB, E) can be extracted from the analysis of the cases (+, — ) and 
(+, +), respectively. 

For the numerical part of this investigation we restrict ourselves to the case N = 1 which 
is the most interesting one in view of applications of the results to binary liquid mixtures. 
The simulations are performed for a spin - | Ising model confined to a film geometry in 
d = 3 dimensions defined by the Hamiltonian 

H! = -J £ s(x,z)s(xV)-iJi]>>(x,i)-# 2 £ S (x,L), (2.3) 

<(X,Z),(X',2')> x x 

where J is the excange coupling constant, < (x, z), (x', z') > denotes a nearest neighbor 
pair of spins and the spins s(x, z) can take the values 1 and —1. The underlying lattice is 
supposed to be simple cubic with L' lattices sites and periodic boundary conditions in the 
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x and y directions. In the z direction the lattice has L < L' sites and the missing bonds 
in the two surface layers at z = 1 and z = L are left open. In order to simulate the model 
at the normal transition Eq. ( |2.3| ) contains two surface terms by which the spins in the two 
surface layers are coupled to surface fields Hi and H 2 , respectively. Infinite surface fields are 
simply realized by fixing all spins in the surface to a fixed value 1 or —1 depending on the 
sign of the surface field. In the model defined by Eq. (|2.3|) the surface exchange coupling 
constant J\ has the fixed value J\ = J. It has been shown by Monte-Carlo simulations of 
spin - | Ising models that the SB multicritical point is characterized by the special value 
J lc ~ 1.50 J of the surface coupling constant J\. Apart from corrections to scaling the 



O surface universality class is represented by the condition J x < J lc which is fulfilled by 
Eq. ( [2. 3D due to J\ = J < J\ c . Therefore only the O surface universality class (Hi = 
or H2 = 0) and the E surface universality class (Hi 7^ or H2 7^ 0) can be studied with 
the above Ising model Hamiltonian. The film geometry underlying Eq. Q2.3|) then allows 
the investigation of the four combinations (0,0), (0,E), (+,+), and (+, — ) of boundary 
conditions by a Monte-Carlo simulation, where the combination (O, E) means (O, +) or, 
equivalently, (O, — ). The principal setup of a Monte-Carlo algorithm for a measurement of 



the Casimir force in lattice models is described in Ref. |33[ to which the reader is referred 
for further details. 

III. LANDAU THEORY 

The presence of a symmetry breaking surface field implies a nonvanishing order parameter 
profile for all r (see Appendix A), which substantially complicates the field theoretic analysis 
of the Casimir effect as compared to the case of symmetry conserving boundary conditions 
discussed in Ref. \$8 \. On the other hand the leading (mean field) contribution to the 



Casimir amplitude can be determined without any detailed knowledge about the functional 
form of the order parameter profile. We briefly illustrate this for the case r = and N = 1, 
i.e., $ = ($,0, ... ,0) in Eqs. (|2.1|) and (|2.2|). In the mean field approximation the order 



parameter profile has the form $(x, z) = (M(z), 0, . . . , 0), where M(z) solves the Euler 
Lagrange equations given by Eqs. ( |A1| ) and (|A2p. Inserting <E» into Eqs. ( |2~T| ) and ( p.2| 



for r = and integrating by parts using Eqs. (|Al|), (|A2|) , and ( |A4] ) H[&] can be evaluated 
without solving the Euler-Lagrange equations for M(z) explicitly. The result is the mean 
field free energy of the film at bulk criticality and is given by 

Km = ^M 2 (0) + |M 2 (L) - hiM(0) - h 2 M(L) 

O O 6 6 

(3.1) 

L 

3 



+ 



l -M'\z») - |M 4 (z ) 



where hi and h 2 denote the first component of hi and Ii2, respectively, and zq is an arbitrary 
reference point < zq < L between the two surfaces of the film. The terms in the first 



line of Eq. (3J.) constitute the surface contribution to the mean field free energy and the 
contribution in the second line of Eq. ( |3.1| ) is the finite size part, where the square bracket 
yields the Casimir force (see below). As a direct implication of Eq. (|A4|) one finds that the 



above expression for the Casimir force does not depend on the reference point Zq. Note that 
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due to t = the bulk contribution to Eq. ( [3.1|) vanishes identically. For r^O 7^[3?j cannot 
be expressed in the same closed form as given by Eq. ( p.l| ) and we therefore resort to the 
zz-component Tj_^ of the stress tensor T^(x, z) in order to find a more general expression 
for the Casimir force. The stress tensor Tjy is given by p|| 



9$ 



Tj Okl 

OXi 



I( V $) 2 + -$ 2 + 4($ 2 ) 2 
2 V ' 2 4! v . 



(3.2) 



d-2 
4(^T 



+ 0(g 3 



d 2 



dx k dxi 



4,.V 2 



where r and g have the same meaning as in Eq. Q2.1|) . The scaling dimension of Tki is given by 
the spatial dimension d. In a film geometry (Tki) is diagonal due to the lateral translational 
invariance of the film. From the conservation property d(Tf-i) /dx k = one then concludes 
that (Tki) does not depend on position and therefore (T±±) can be directly identified with 
the Casimir force per unit area. Note that the evaluation of (Tjy) according to Eq. ( |3.2| ) for 
Xk = xi = z and for r = within the mean field approximation $(x, z) = (M(z), 0, . . . , 0) 
for the order parameter yields the square bracket in Eq. ( p.l|) . 

We now turn to the mean field analysis of the Casimir force as a function of the reduced 
temperature r, where we first restrict ourselves to the case N — 1 (Ising universality class). 
In view of later applications of the results to binary liquid mixtures near the critical demix- 
ing transition this is the most relevant case. For the mean field analysis alone it would not 
be neccessary to determine the full order parameter profiles. However, in order to perform 
the fluctuation expansion (see Sec. IV and Appendix B) precise knowledge about the pro- 
files on the mean field level is indispensable. Details of the calculation are summarized in 
Appendix A. In the course of the calculations for the order parameter profiles one obtains 
the corresponding expressions for the Casimir forces as byproducts which will be discussed 
in the following paragraph. 

As in Appendix A we write the mean field contribution (Tj_j_) to the Casimir force in 
the form (Tj_j_)o = (Q/g)t±± and we only consider t±± in the following for simplicity. From 
the general theory of critical finite size scaling |7]|J we expect tj_± to take the scaling form 

(3.3) 



t±± = L~ d F{y) , y = tL 1 '^ 



where d = 4 and v — 1/2 within mean field theory. Note that right at the upper critical 
dimension d = d c = 4 the prefactor 6/ g of the Casimir force generates logarithmic finite-size 
corrections due to the fact that the renormalized counterpart u of the coupling constant g 
vanishes according to u(l) ~ 1/lnZ for / — > oo at the renormalization group fixed point j45 . 
However, logarithmic corrections to scaling in d — 4 will be disregarded here so that from 
the point of view of mean field theory the above prefactor is treated as a constant. 

For the case of (+, +) boundary conditions the scaling function F +i+ (y) can be read off 
from Eqs. QAUD and ([AT]]). The result is 

(2K) 2 (2k 2 
-(2K) 2 (k 2 



F +>+ (y) = -(2K)*k 2 (l-e) 



4 7.2 



F + , + (y) = 



y 
y 



- 1); 
+ D, 



(3.4) 



where K = K(k) is the complete elliptic integral of the first kind and < k < 1. The y 
dependence of F +t+ according to Eq. (|3.4j ) is given in the parametric form y = y(k), where 
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y(k) is a monotonic function of k so that the inverse k = k(y) exists and constitutes the y 
dependence of F +t+ in a unique way. As can be seen from Eq. (|3.4|) the parameterizations of 
F +i+ and y for y > —tt 2 and y < — tt 2 are different. The reason for this is purely technical in 
the sense that negative values for k 2 are avoided (see Appendix A for details). There is no 
singularity of F +t+ (y) at the point y = —tt 2 (k = 0). In fact, F + + (y) is analytic for all values 
of y, because the critical point of the film (T C (L), h c (L)) is located off the temperature axis 
at a finite critical bulk field h = h c (L) ~ L~ A / U , where A is the gap exponent PJITJ. Within 



mean field theory one has A = 3/2 so that A/V = 3. For (+, — ) boundary conditions the 
corresponding result for the scaling function F +r _(y) can be read off from Eqs. ( |A19 ) and 
( ggg ). One finds 

F + ^y) = {2K)\l-k 2 ) 2 , y = 2{2K) 2 {k 2 + l) ] 

F + ^{y) = {2Kf , y = -2(2K) 2 (2k 2 -l), (3.5) 



where a parameterization analogous to the one in Eq. ( |3.4|) has been used. The scaling 
function F +i _(y) is also analytic for all values of y, although the critical point of the film 
in the case of opposing surface fields is located on the temperature axis and is associated 
with the interface derealization transition [|TBJ. However, due to the limit hi = — h 2 — > oo 
performed here this critical point has been formally shifted to y c = — oo so that it is no longer 
visible as a singularity in F + _(y). Corresponding results for (SB, +) and (O, +) boundary 
conditions can be constructed from Eqs. ( |3.4| ) and (|3.5| ) using the simple transformation 
L — > 2L (see Appendix A). One obtains 

F SB ,+{V) = T6 F +,+M , Fo,+(v) = Te F +A±y)- (3-6) 

The scaling functions obtained so far still contain a bulk contribution which corresponds to a 
bulk pressure given by t±± t buik = —rml—mf. For r > one has t±^ tbu [ k = {m b = 0) and for 
r < one has t±± jbu i k = r 2 /4 (m h = ^J—r/2). The bulk contribution F bu i k (y) to the scaling 
functions given by Eqs. (|3.4j), (p.5|), and ( |3.6| ) then has the simple form F buUt (y) = 9(—y) y 2 /4: 
which contains the usual mean field singularity of the bulk free energy at r = 0. In order to 
express the finite-size contribution to the Casimir force (Tj_j_)q in the scaling form we define 
the scaling functions 

faM = F aM ~ F bulk{y) (3.7) 

which are displayed in Fig. [l| for (a, b) = (+, +) and (+,—). Their shapes resemble those 



of the corresponding scaling functions for the Ising model confined to a strip in d = 2 [|30 
The asymptotic behavior of the scaling functions for y — > ±oo is governed by an exponential 
decay according to 

f+Av -> oo) - - 1 % 2 exp(- v /y), 
f+,+{v -oo) ~ -16y 2 exp(-^2y), 

f+-(y o°) - exp(— v^), 
.2 



/+,_(</ - -00) ~ 16y 2 exp(-^/-y/2). (3.8) 

The scaling functions take quite sizable values over a surprisingly broad range of the scaling 
argument y. This may serve as a first indication that the Casimir forces provide a strong 
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modification of the usual dispersion forces in a parallel plate geometry at the extraordinary 
transition. However, in order to estimate the absolute strength of the Casimir forces in 
binary liquid mixtures close to their critical demixing transition a renormalization group 
analysis of f +!+ (y) and f+-(y) is required (see Sec. IV). 
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Scaling functions /+ 



(y) (solid line) and f+ t -(y) (dashed line) according 
(3.5), and (|3.7|) . The y range influenced by the bulk critical point 
is very broad and the asymptotic decay for y — > ±oo is dominated by an 
exponential [see Eq. (|3.S|)1. Note that /+ t +(y) and f+-(y) take their extreme values 
at y ~ 10 and y ~ —25, respectively. 



If the order parameter has N > 1 components, the case of parallel surface fields is 
already covered by the above analysis of the (+, +) boundary conditions for N = 1, because 
in this case the order parameter only has one nonzero component parallel to the surface 
fields. For antiparallel surface fields, however, this is not as obvious, because the order 
parameter has the additional freedom to rotate across the film by a position dependent 
angle f(z). We illustrate this for the case N > 2 with h x = [hi, 0, . . . , 0) and h 2 = 
hi (cos a, sin a, 0, . . . , 0) in the limit hi — > oo and for r = 0. A similar situation has been 
been discussed in Ref. |39[] for r ^ 0. If the order parameter profile is written in the form 
Wl(z) = <JV2/g m(z) {cos(p{z), sin</?(z), 0, . . . , 0) one finds the Euler-Lagrange equations 
[see Eq. (fa) and Ref. |1] 



(p'(z)m 2 (z 

m"(z) 



0, 

(f' 2 (z)m(z) + 2m 3 (z). 



(3.9) 



As boundary conditions for <p(z) we choose <£>(0) = and (f(L) = a, because the order 
parameter should be parallel to h x and h 2 , respectively, at the surfaces. The amplitude 
function m(z) is positive and its qualitative behavior resembles that of the profile m +t+ (z) 
[see Eqs. ( |A13|) and (|A 1 5|) ] . From Eq. ( |3.2|) we then find for the Casimir force 



S 



(7± ± )o = (Q/g) t ±± = (6/g) ^(L/2)m 2 (L/2) - m\L/2) 



(3.10) 



Note that Eq. ( |3.10| ) allows a sign change of the Casimir force as a function of the angle 
a enclosed by the surface fields. Following Appendix A [see Eqs. (|A4|) and (|S5|)] the first 
integral of Eq. (|3T9|) is given by 



m' 2 (z) 



c/m 2 (z), 

-c 2 /m 2 (z) +m 4 (z) +tj_i_, 



where 



t±± = c 2 /m 2 {L/2) - m\L/2) 



and c is a constant such that 

tp(L) for (p(z) = c dz' /m 2 (z'). 



a 



(3.11) 



(3.12) 



(3.13) 



Just as for Eq. ( |A6| ) we apply the substitution P(z) = m 2 (z) and eliminate c using Eq. 
( i.1'2 ). All the information needed to calculate the Casimir force, i.e., t±± as a function of 
a is now contained in Eq. ( |3.13| ) and 



P' 2 (z) = 4 P 3 (z) - P\L/2) + t ±± (P(z) - P(L/2)) 



(3.14) 



which shows that P(z) = p(z; <7 2; <?3) is a Weierstrass elliptic function, where the invariants 



#2 and #3 can be read off from Eq. ( 3.14 ). As we are focussing on the limit of infinite surface 
fields the film thickness L is one of the basic periods of P{z) [see Eqs. ( |A1 1| ) and ( |A18 )1, 
and therefore P{z) has double poles at z = and z = L. Using Eq. ( |3.14|) we can rewrite 
Eq. ( |3.13| ) and find a representation for P(L/2) by performing a separation of variables in 
Eq. (|3.14|) . Writing t±± in the scaling form [see Eq. (|3~3|)1 



t±± = L 4 g(a) 



(3.15) 



and using the abbreviation p = LJP(L/2) one finds 



V 



x 3 - 1 + (x - l)p 4 g(a) 



-1/2 



dx. 



(3.16) 



a — Jl +p~ 4 g(a) / x 1 x 3 - 1 + (x - l)p 4 g(a) 



-1/2 



dx. 



The solution of Eq. ( |3.16|) is shown in Fig. ^|. The Casimir force (i.e., g{a)) grows monoton- 
ically from a = to a = ir at fixed L and vanishes for the angle a = tt/3 which can also be 
derived directly from Eq. ( |3.16|) by setting g (a) = 0. Furthermore it should be noted that 
according to Eqs. ( |3.4p and ( |3.7| ) one has g(0) = F + + (0) = / +;+ (0) and according to Eqs. 
( |3.5| ) and ( p.7| ) one also has g(ii) = F +i _(0) = / +i _(0). The function g(a) therefore smoothly 
interpolates between (+, +) and (+, — ) boundary conditions giving the same result for the 
Casimir force as the Ising universality class (N = 1) in these two cases. 
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FIG. 2. Amplitude function g(a) of the Casimir force according to Eqs. ( |3.15| ) 
and ( 3.16| ). g(a) smoothly interpolates between g(0) = /+ i+ (0) and g(Tr) = / +) _(0) 
(see main text). The amplitude vanishes at a = tt/3. 



We close this section with a short discussion of the analytic solution of Eqs. ( |3.14p and 
( p,16| ). Following the derivation described in Appendix A and using Eq. ( |3.13|) the profile 

rn(z) = \Jp(z), the amplitude function g(a), and the angle a can be parameterized in terms 
of the modulus k of the Jacobian elliptic functions. One finds 



m(z) 



2K 



1 



g(a) = --(2Ky 



dn 2 {(-k) 2k 2 -1 
sn 2 (C; k) + 3 

(2k 2 - l) 2 



1/2 



1 + 



(3.17) 



a = 2 



■\-2k 2 )(2-k 2 
3(1+ A; 2 ) 



1/2 



n(l/3 + k 2 /3, k) - K 



and 



2K 



m(z) 

g(a) = [2Kf 



cn 2 ((;k) 



2{2k 2 



.1/2 



L [sn 2 (C;A;)dn 2 (C;A;) 

A(2k 2 - l) 2 



1 - 



(3.18) 



a = JQ{l-2k 2 ) 



1 - 2k 2 + 3 



[n(a, k) - K] 
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[U(b, k) - K] 



l-2k 2 

where ( = (2K/L)z, < k 2 < 1/2, and the parameters a and b are given by 



a 
b 



1 + 4k 2 + J9 - 8(2k 2 - 1) 



1 + Ak 2 - J9-8(2k 2 -l) 



(3.19) 



Furthermore K = K(k) and H(x, k) for x = a,b denote the complete elliptic integrals of the 
first and the third kind, respectively. The angle a traverses the interval [0, tt] by decreasing 
k 2 from k 2 = 1/2 to k 2 = in Eq. (|3~T7|) , changing to Eqs. ( P^ ) and ( ET19D at A; 2 = 
and increasing /c 2 back to k 2 = 1/2. The special point k 2 = has no particular physical 
significance, it only marks a singular point in the above parametric representation of g(a). 
From Eq. (|3.18|) we identify k 2 = (2 — v / 3)/4 as the parameter value, where g(a) = or, 
equivalently, a = ir/3. 




FIG. 3. Amplitude m(z) (solid line) and phase (p(z) (dashed line) of a two 
component order parameter for surface fields at an angle of a = 0.98-7T [k 2 ~ 0.499, 
sec Eqs. ( |3l3l) and (|3Tl8|) and main text]. m(z) has been normalized to m(L /10) 
so that m(z) and <p(z) can be plotted on the same scale. 



Setting k 2 = 1/2 in Eq. (|3T7D yields m(z) = m +j+ (z) [see Eq. (I&T3D for r = 0], 
= 0) = = 0) [see Eqs. (Q and Q], and <p(z) = [see Eq. (&Wj)}. This means 

that the order parameter profile is given by M(^) = yl2/g (m +t+ (z), 0, . . . , 0) as anticipated 
from the case N — 1 for (+, +) boundary conditions. In the limit k 2 1/2 Eqs. (|3.18| ) and 
( gl9| ) yield m(z) = |m+,_(z)| [see Eq. (|A20|) for r = 0], g(a = vr) = = 0) [see Eqs. 



( P-5| ) and (|3.7|) 1, whereas v 9 (- 2 ) here is given by the step function <p(z) = n 9{z/ L — 1/2). The 
order parameter profile is then given by M(z) = \Jl2/ g(m + ^(z), 0, . . . , 0) which shows that 
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also for antiparallel surface fields mean field theory for an A-component order parameter is 
already captured by the case N = 1. We illustrate this remarkable behavior of M(z) for 
a/ir = 0.98, i.e., a situation close to antiparallel surface fields. The corresponding modulus 
k [see Eqs. ( [3.18| ) and (|3.19| )1 is given by k 2 ~ 0.499. The phase <p(z) and the amplitude 
m(z) of the order parameter are shown in Fig. |3|. The order parameter rotates by almost 
the full amount a in a narrow interval around z = L/2, where m(z) is smallest. In the limit 
a —>■ ix this interval shrinks to the point z = L/2, where m(z) vanishes and (p(z) becomes 
discontinuous. 

Although the Casimir force is governed by universal scaling functions IPSlfl it is not 
possible to estimate their absolute magnitude within Landau (mean field) theory. The 
reason is that for the boundary conditions considered here these scaling functions contain a 
common prefactor which depends on the bare coupling constant g and therefore has a value 
inaccessible by pure mean field arguments. In order to at least partly fill this gap we now 
turn to the field theoretic analysis of the Casimir force at bulk criticality. 



IV. CASIMIR AMPLITUDES 



At the bulk critical temperature T = T c (r = 0) the Casimir forces in a film are governed 
by the universal Casimir amplitudes which explicitly depend on the two surface uni- 
versality classes combined in the film. For N > 1 order parameter components the Casimir 
amplitudes may also depend on continuously varying parameters, as demonstrated above 
for the case N = 2 with tilted surface fields. For T ^ T c (r ^ 0) these amplitudes have 
to be replaced by universal scaling functions Q a ,b(y) of a suitably chosen scaling argument y 
28| , which will not be considered in this section. 

The Casimir amplitude is defined as the finite-size amplitude of the free energy of a film 
at bulk criticality [P,EH] . Translating this definition to the force one finds 



(Tlx) = (d-l)A a , b L 



-d 



(4.1) 



in d dimensions and for r = 0, where f(T c , L) = TiT = T c , L)/(A ksT^) is the critical part 
of the free energy per unit area A of the plates. Following Ref. ksT c is used as the 
natural energy unit for the free energy, where k% is the Boltzmann constant. As a first step 
beyond Landau theory the contribution of Gaussian fluctuations to the Casimir force, i.e., 
to the amplitudes A a & will be investigated here. We introduce the fluctuation part <fr of the 

order parameter <fr by <3> = $ + where <fr = M(z) = \Jl2/g (m(z), 0, . . . , 0) is the mean 
field order parameter profile discussed in the preceding section and Appendix A. Inserting 
the above decomposition of $ into Eqs. (|2.1|) and ( pT2|) for r = and C\ = c 2 = and 
keeping only the quadratic terms in $ = ((f) i, . . . , (f>^) we obtain 



=H[®] + ^Jd d - 1 x£dz{(V<f>) 2 + [6m 2 (z)(Pl + 2m 2 (z) ($ + ... + &)]} 



+ o 



:$ 2 ) 2 



(4.2) 



where hi = ±h 2 = (hi, 0, 
contribution to Eq. 



0) in the limit hi —>■ oo is implicitly assumed. The mean field 
has already been discussed in Eq. ([3.1|). Following Eq. (IO) 
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we decompose the Casimir force into the mean field part, a Gaussian part, and higher order 
corrections according to 



(T±±) = (T ±± ) + (T ±± ) l + 0(g) = (Q/g)t ±± + (T ±± )x + 0(g). 



(4.3) 



In order to determine (Tj_j_)i from Eq. ( |3.2| ) one also needs the cubic terms in Eq. ( |4.2| ) 
and we will therefore not follow this approach any further. It is much more convenient 
to determine (Tj_j_)i from the Gaussian contribution to the free energy by taking its first 
derivative with respect to the film thickness L p8 |. Following Ref. |28| this can be done 
most easily in a spectral representation of the Gaussian Hamiltonian given by Eq. (|4.2| ). For 
the evaluation of only the eigenvalue spectrum is needed. According to Eq. (|4.2|) the 

spectrum consists of a longitudinal part characterizing the eigenmodes of the longitudinal 
fluctuations 0i of the order parameter and a transverse part e^ 1 ^ which is the same for each 
of the N — 1 transverse components (02; • • • > 4>n) of the order parameter fluctuations. The 
spectra eW and effi are determined in Appendix B. Once the eigenvalues are given one can 
employ the dimensional regularization scheme and according to Ref. we find 



f(T c ,L) = nm 



r [(3 - d)/2] 

2 <f-i 7r (d-i)/2( d _ 



oo 
n=3 



(d-l)/2 



+ (N 



oo 

E 

n=2 



(*>) 



(d-l)/2 



(4.4) 



for the critical part free energy within the Gaussian approximation in d dimensions. The 
L dependence of the Gaussian contribution to f(T c , L) is completely determined by the L 
dependence of the eigenvalues. From simple dimensional analysis one has ~ L~ 2 so that 
defjdL = -(2/L)ef for i = 1, 2. From Eqs. (O) and (O) we find 



(T±±) 



6 
-1 

9 



-t±± 



r [(3 - d)/2] 

2 d-l 7r (d-l)/2 L 



E(4 2) ) (d " 1)/2 + (^-DE(^ 



(d-l)/2 



.n=3 



+ 0(g) (4.5) 



for the Casimir force in the Gaussian approximation. The mode sums in Eq. ( |4.5|) diverge for 
d = 4 — e and we therefore employ the dimensional regularization scheme. Furthermore, the 
above sums yield an UV singularity in the typical form 1/e which must be treated analytically 
in order to facilitate the renomalization of Eq. (|4.5| ). Both objectives can be achieved with 
the asymptotic expansions of the eigenvalues effl and for large mode numbers n which 
are given by Eqs. ( [B 1 3|) and ( B15| ). The regularization of the mode sums and the analytical 
treatment of the 1/e pole is summarized in Appendix C. Using the results from Appendix B 
and C we now investigate the different boundary conditions separately, where the mean-field 
results given by Eqs. (|3.4p, (|3.5|), and ( |3.6| ) are only needed for k 2 = 1/2 (r = 0). 

For the renormalization of the Casimir force given by Eq. ([4.5|) we use the conventions 



of Ref. [ 28 1 and define the renormalized coupling constant u by 



g = 2 7r 1 n Z u 



N + 8u 
3 1 



0(u 2 



(4.6) 



where /i is an arbitrary momentum scale and d = 4 — e in the following. The infrared stable 
fixed point value u*(e) of the renormalized coupling constant u is given by 



3 9(3iV + 14) 2 . , 3 . 



(4.7) 
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For later reference we also quote the 3-loop estimate [f46 



e 17 2 
3 + + 



709 
17496 



- C(3) ) + 0(e* 



(4.8) 



of the fixed point value u* for iV = 1, where ((3) ~ 1.20206 is a special value of the Riemann 
zeta function. In order to improve the predictive quality of a low-order e-expansion for e = 1 
(d = 3) in a simple way one may try to include exact results for the quantity in question 
in d = 2 in the sprit of a Pade approximant in the variable e. This can be applied rather 
successfully to the Casimir amplitude A per thus improving the agreement between the field- 
theoretic prediction |2£| and the Monte-Carlo estimate |33| in d — 3. We will therefore follow 
the same procedure here, where the case of (SB, +) boundary conditions must be excluded, 
because the SB multicritical point does not exist in d = 2. 

The renormalized expression for (T±±) for (+, +) boundary conditions can be obtained 
by inserting the mean field result given by Eq. Q3.4| ) for y = and the regularized mode sum 
given by Eq. ( |C5|) into Eq. ( |4.5|) and by applying the renormalization prescription given by 
Eq. (t4.6|) . After expanding all d dependent quantities to first order in e [see Eq. (|C12j )1 the 
1/e pole coming from Eq. (|C5| ) is cancelled, i.e., the UV singularity has been consistently 
removed from the theory. The Casimir force then follows by evaluating the resulting renor- 
malized expression for (T±±) at the renormalization group fixed point u = u*(e) given by Eq. 
( ]4.7|) . The e-expansion of the universal Casimir amplitude A + + , which characterizes the 
strength of the Casimir force in a critical film with parallel surface fields, is finally obtained 
by applying the definition of A aj b given by Eq. ( }4.1|) . The algebraic manipulations involved 
here starting from Eqs. ( |4.5| ), flC5]), and ( f4.6| ) are absolutely elementary, so that we only 
quote the final result 



A, 



7T d / 2 T{d/2) /KV 
2u*(e) V 7r ) 



0.6853 + e 0.1242 + 0(e 2 

N+8 N+8 K 



(4.9) 



where part of the e-expansion has been resummed consistently to first order in e using Eq. 
( pi3| ). For (+, — ) boundary conditions the same procedure can be applied using Eq. ( p.5| ) 
for y = and Eqs. fl4.5|) and ( |C6|) . One finds 



A_ 



7T 



d / 2 T(d/2) (K\ 



(4.10) 



From Eq. ( |3.6| ) for y = and Eqs. 

n d / 2 T(d/2) /K s 4 



9f /V — 1 

— — 0.2822 + e 0.4066 + 0(e 2 ) 

N + 8 N + 8 K ' 

and ( |C9D one has for (SB, +) boundary conditions 



ISB, 



32u*(e) 



71 



9e N — 1 

1 + t^—: 1-7141 + e ^ 2.8448 + 0(e 2 ) 



N 



N + 8 



(4.11) 



From Eq. ( |3.6| ) for y — and Eqs. ( |4.5| ) and (|C11| ) one finally has for (O, +) boundary 
conditions 



7T d / 2 T{d/2) /K 
8u*(e) V 7T 



•1 r 



1 + 



9£ 



A^ + 8 



0.1988 + e 



N - 1 

ivTs 



0.2289 + 0{e 2 ) 



(4.12) 



Note that u*(e) in the above expressions is given by Eq. ( f4.7|) . It is remarkable that the 
coefficients of the Gaussian contribution to ^sb,+ given by Eq. (|4.11| ) are much bigger than 
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the corresponding coefficients in Eqs. ( |4.9|) , (|4.10|) , and ( [4.12|) . This may be due to the fact 
that the order parameter near a SB surface is much more susceptible to fluctuations than 
near O or E surfaces. 

In the Ising universality class (N = 1) in d = 2 three of the above Casimir amplitudes 
are known exactly from conformal field theory |29[. They are given by 



A_ 



n 23 

, A , = — 7T 

48 + ' 48 



71 



(4.13) 



The construction of a Pade approximant from Eqs. (|4.9| ), G4.10|) , and 12|) for N = 1 
which extrapolates to the amplitudes given by Eq. ( 4.13|) for e = 2 is arbitrary to a certain 
degree. If one uses Eq. fl4~8| ) instead of Eq. (|4.71 ) for iV = 1 and introduces an additional 
^-contribution to the square bracket of Eqs. (f4.9|) , ( |4.10| ), and ( |4.12| ) such that Eq. (|4.13| ) 
is reproduced for e = 2 one finds the interpolation formulas 



A + _ 



7i d ' 2 T(d/2) (K_ 
2u\{e) \~ 

7i d ' 2 T(d/2) (K 



7T 



A 



o, 



7T d / 2 T{d/2) fK_ 



1 - 0.6853 £ + 0.1275 e x 



l- 0.2822 e + 0.0914 t 



1 + 0.1988 e- 0.0707 e' 



(4.14) 



Numerical estimates of the Casimir amplitudes in d 
formulas are summarized in Table I. 



3 obtained from the above analytical 



Table I. Casimir amplitudes for the Ising universality class in d = 3. The values labelled 
by e = 1 are obtained by evaluating Eqs. (|4.17[) , (p|), |TlCl ), (|4TT1) , and (|1|) for iV = 
1 and e = 1. The values labelled d = 3 are obtained from Eqs. ( [4.18 ) and ( f4.14| ) for d = 3 
(e = 1). The Monte-Carlo estimates obtained from the serial version of the algorithm 
presented in Ref. |33| are labelled by 'MC (see also main text). Statistical errors (one 
standard deviation) are in the last two digits as indicated inside the parenthesis. The 
last line shows Migdal-Kadanoff estimates taken from Ref. ||3z 





Aper 


A ,o 






Ass,+ 


Ao,+ 


e = 1 


-0.1116 


-0.0139 


-0.173 


1.58 


-0.093 


0.165 


d = 3 


-0.1315 


-0.0164 


-0.326 


2.39 




0.208 


MC 


-0.1526(10) 


-0.0114(20) 


-0.345(16) 


2.450(32) 




0.1873(70) 


Ref. 


I 


-0.015 





0.279 


0.017 


0.051 



A, 



and A_i_ _ can be 



In d — 3 and for N — 1 the Casimir amplitudes Ao,o, Ao, 
measured by a Monte-Carlo simulation of the Ising model defined by Eq. ( |2.3| ). The algo- 
rithm and its special adaptation to the measurement of the Casimir amplitude is presented 
in Ref. [33] in detail. We therefore only briefly describe the differences between the imple- 
mentations used here and in Ref. |33|. The present implementation of the algorithm utilizes 



a serial hybrid update scheme which consists of a Metropolis update sweep of the whole 
lattice followed by a Wolff update. The length of the equilibration and the measurement 
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period used here correspond to those in Ref . |33| . The slab geometry contains L' 2 x L lattice 
sites, where L'/L must be chosen as large as possible in order to approximate the infinite slab 
geometry. In practice already L' = 4L turns out to be sufficient, i.e., the results obtained 
for this choice agree with those for V = QL within a fraction of one standard deviation. The 
thickness L of the slab has been varied between L — 12 and L = 32 layers. As in Ref. |!3| we 
use the multiple histogram technique |47] , where the number of histograms taken has been 
increased from 25 to 31 for L > 24 on order to guarantee sufficient overlap between adjacent 
histograms |33] . The simulations were run on DEC Alpha workstations at the University of 
Wuppertal and the total amount of CPU time used is equvalent to about one year of CPU 
time on a DEC 3000 workstation. 

The serial implementation of the algorithm has been tested for the Casimir amplitude 
A per with V = 4L and L' = 6L for L = 20 and L = 24. The estimates for A per obtained 
with these four lattice sizes agree within their statistical error and give the final estimate 



A 



per 



-0.1526 ±0.0010, 



(4.15) 



which is in perfect greement with the estimate obtained from the parallel algorithm ||33|| . 
The amplitude Ao,o has been measured for the same lattices sizes and for Ao,+ additional 
simulations were performed with L = 28 and L' = 4L. All individual measurements agree 
within their statistical error and the final estimates are shown in Table I. For (+, +) and 
(+, — ) boundary conditions, however, the situation is different. For A + + measurements 
have been made for 12 < L < 32 and V = 4L, the individual estimates are displayed in 
function of L. 
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FIG. 4. Monte-Carlo estimates of the Casimir amplitude A +)+ as a function of 
the number of layers L in a L' 2 x L slab for V = 4L. The size of the error bars 
represents one standard deviation. The data point at L = 32 is taken as the final 
estimate (see Table I). 



The estimates show a clear systematic dependence on L and apparently even for L = 32 
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layers the asymptotic regime has not yet been reached. The last three data points fall onto a 
straight line within their error bars so that the data cannot be extrapolated to an asymptotic 
value. As the current Monte-Carlo estimate for A + + we therefore take the measurement 
for the biggest system (L = 32, L' = AL) (see Table I). The situation for A +j _ is similar. 
The individual measurements are shown in Fig. ^ for 12 < L < 28 (L' = AL). Again, the 
asymptotic regime has not been reached for the biggest system, but this time it is possible 
to estimate the asymptotic value for A +i _ by a least square fit of the function 



A e Jf_{L) = A +i _ + Dexp{-KL) 

to the data for L > 16 using A + _, D, and k as fit parameters. 
2.40 



(4.16) 




FIG. 5. Monte-Carlo estimates of the Casimir amplitude A +j _ as a function of 
the number of layers L in a L' 2 x L slab for V = AL. The size of the error bars 
represents one standard deviation. The dashed line shows a fit of Eq. ( |4.16 ) to the 
data for L > 16 giving the estimate of A +j _ displayed in Table I. The absolute size 
of the errorbars is about twice that in Fig. 0. 



The exponential L dependence of A+T(L) in Eq. ( ^4.16 ) is motivated by the short- 
ranged nature of the interaction in Eq. ( |2.3j ). The error of the amplitude A +i _ is estimated 
by taking the maximal error of the individual measurements involved in the fit. All estimates 
obtained from Eqs. (^9[), ( |4.iq ), and ( f4.12| ) for N = 1, from Eq. ( |4.14| ), and our Monte- 
Carlo estimates are summarized in Table I. For completeness we also display estimates for 
A„ er and Ao,o obtained from the partially resummed e-expansions [ETI] 



A 



per 



-N 



r(d/2)C(d) 



7T 



d/2 



5N + 2 
AN + 8 



e + 0{e 2 



A 



o,o 



2~ d A 



per 



(4.17) 



for N = 1 and from the Pade approximants ||2T| . |33 ] 
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\er — ^[jt~ " y- ~ 4 7 _ 1 ) ' ^o,o — 2 d A per , (4-18) 



which reproduce the exact results |29 



A p er = -^ , A 0j0 = -^ (4.19) 

in d = 2. For comparison we also reproduce Migdal-Kadanoff estimates for the Casimir 
amplitudes in d = 3 from Ref. ||32|| . The agreement between the Pade approximants and 
the Monte-Carlo estimates is quite satisfactory, except for Ao,o which seems to be closer 
to the partially resummed e-expansion and the Migdal-Kadanoff estimate. However, the 
amplitude is rather small and therefore the relative statistical error of the Monte-Carlo 
estimate, which is one standard deviation, is very large (20%, see Table I). In view of Fig. 
f| the Monte-Carlo estimate for A + + given in Table I constitutes only an upper bound for 
the true amplitude and must therefore also be handled with caution. The fit procedure used 
to extract A +i _ from the data shown in Fig. |5] is also susceptible to systematic errors to a 
certain extent. However, compared to the parameters D and k in Eq. ( [4.16|) the resulting 



estimate for A + _ is quite robust with respect to, e.g., changes in the number of data points 
included in the fit. The obtained variation of A +j _ is in the same order of magnitude as 
the statistical error given in Table I. With regard to their reliability the analytical and the 
Monte-Carlo estimates of A +j+ , A +i _, and Ao,+ seem to be a substantial improvement over 
the Migdal-Kadanoff results. 

V. EXPERIMENTAL IMPLICATIONS 

A typical experimental setting, within which the film geometry considered here is of 
particular interest, is provided by wetting experiments performed on plane and chemically 
homogeneous substrates P, PH[^H| . The equilibrium thickness L of the wetting layer is de- 
termined by the minimum of the effective interface potential j3|. It is given by the grand 
canonical free energy of a liquid layer of a prescribed thickness I, which is in contact with 
the substrate on one side and with the bulk vapor phase on the other side. In the limit of 



large interfacial areas A the effective interface potential can be written in the form |3,48,49 



lim 9£Efll = u (i) = i [ Pl (T)/p v (T)-l]p {T) 5p + a sl (T) + a lv (T) + 8u(T,l), (5.1) 

where pi(T) and p v (T) are the liquid and the vapor density, respectively and Po(T) denotes 
the liquid- vapor coexistence line in a p, T phase diagram. The quantity dp in Eq. ( |5.1| ) is a 
dimensionless measure of the undersaturation of the vapor, i.e., bp > indicates that in the 
bulk the vapor phase is thermodynamically stable. The substrate-liquid and liquid-vapor 
interfacial tensions a s i{T) and cri v (T) do not depend on / and Sco(T, I) contains the dispersion 
(van der Waals) forces and the critical Casimir forces in the liquid layer. For a binary liquid 
mixture as the wetting agent the critical point of interest is the critical end point of the line 
of critical demixing transitions on the liquid-vapor coexistence surface (see Fig. 1 in Ref. 
||49||). In order to discuss the effect of criticality on the equilibrium thickness L of the wetting 
layer |10| , ^8| we assume in the following that the critical temperature T cep associated with 
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this critical end point is located above the wetting temperature T w so that the condition 
T w < T ~ T cep guarantees a macroscopic wetting layer of a critical binary liquid mixture. 



For large values of I the van der Waals contribution to 6lu(T, I) has the asymptotic form |50 



_ / W(T) I' 2 + 0(l~ 3 ) (nonretarded) 
OcW^i , <j - | ^ (T) z _ 3 + ^_ 4 ) ( retarded ). V-V 

The explicit temperature dependence of the Hamaker constant W(T) and its retarded 
counterpart W r (T) is quite weak and can be disregarded in the critical regime around 
T cep . According to Eq. ( f5.1|) one has with Suj(T,1) = 5u v dw(T,l) taken from Eq. ( |5.2|) 
L(5p) oc (5p)~ l l 3 in the nonretarded case and L(5p) oc (<5p) -1 / 4 in the retarded case. Pro- 
vided, the the wetting layer becomes thick enough, one observes a crossover from the former 
to the latter power law for 5p — > in a wetting experiment, because the van der Waals 
forces become retarded as L increases [j50|1 . At the critical end point 5u is modified by the 



long-ranged Casimir forces according to 

5uj(T cep , I) = 5uj vdW (T cep} I) + fc B T cep A a)6 r (rf_1) (5.3) 

in d dimensions, where k B is the Boltzmann constant and A a ^ is the Casimir amplitude 
for boundary conditions of type (a, b) as discussed in the preceding section. If the van der 
Waals forces are not retarded one can combine Eqs. (|5.2| ) and (|5.3| ) in d = 3 by defining the 
effective Hamaker constant HSf 



W eff = W + k B T cep A a:b , (5.4) 

where the temperature dependence of W has been disregarded. The effective Hamaker 
constant W e ff replaces W in the effective interface potential given by Eq. (|5.1|) and thus 
determines the equilibrium thickness L of the wetting layer for fixed undersaturation Sp. 
The ratio R(Sp) of the wetting layer thickness L cep (5p) at the critical end point and the 
thickness L(5p) of the wetting layer outside the critical regime is then determined by the 
ratio W e ff/W H|. One obtains 



R(Sp) = L cep (6p)/L(6p) = (W eff /Wy/ 3 = (1 + k B T cep A a , b /wy/ 3 , (5.5) 

which is independent of the undersaturation 5p to leading order in 5p (see Ref. [Q for de- 
tails). If both the liquid-substrate and the liquid- vapor interface prefer the same component 
of the binary liquid mixture one has (a, b) = (+, +) and Eq. ( |5.5| ) predicts a thinning of 
the wetting layer, because A +i+ < (see Table I). In the opposite case (a, b) = (+, — ) 
applies and Eq. (|5.5|) predicts an increase in the wetting layer thickness due to A +j _ > 0. 



An experimental realization for the latter case is provided by a methanol-hexane mixture 



on Si - Si02 wafers as substrates |[T]. The mixture wets the wafers at a temperature be- 
low T cep ~ 300-fT, where the methanol concentration is enhanced near the substrate and 
the hexane concentration is enhanced near the liquid-vapor interface providing a realiza- 
tion of the (+, — ) boundary condition. The Hamaker constant for this system is given by 
W ~ 9 x 10~ 15 erg |51] and with A + _ ~ 2.4 taken from Table I one obtains R(6p — >• 0) ~ 2.3 
from Eq. ( |5.5| ). The corresponding value of R for 4 He on Ne substrates at the lower A-point 



is R ~ 0.995 [48]. The explanation for this drastic difference is twofold. First, there is the 
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combined effect of the Hamaker constant W and the relevant energy scale given by /cbT c . 
For methanol-hexane on Si - Si0 2 one has T c = T cep ~ 300K so that W / '(/csT cep ) ~ 0.2, 
whereas for 4 He on Ne one has T C = T\ = 2.17K which implies W/ (ksT\) ~ 2 ||48|| . Second, 
the relavant Casimir amplitude is A +> _ ~ 2.4 for methanol hexane and Ao,o — —0.022 
for 4 He |[48|| . In the ratio (W e ff — W)/W [see Eqs. ( |5.4j ) and ( |5.5| )1 one therefore has one 



factor ~ 10 in favor of methanol-hexane coming from WjiksTc) and a second factor ~ 100 
in favor of methanol-hexane from the Casimir amplitude which combine to the observed 
drastic quantitative difference in R(5p). 

For 5p — > the equilibrium thickness L(5p) of the wetting layer increases so that the van 



der Waals forces may become retarded [see Eq. (|5.2j) l . In the retarded regime the critical 
contribution to 5u>(T cep , I) becomes the leading term in Eq. ( |5.3|) for d = 3 and therefore 
R(Sp) defined by Eq. ( |5.5| ) diverges for Sp — > according to 

R(Sp _ 0) = (U^A V3 ^ (Pl) ^ {5p) -Vu (5 . 6) 

V Pl-Pv J \Pl-PvJ \PvJ 

For (+, +) boundary conditions one has A +i+ < and in this case retardation of the van der 
Waals forces leads to a finite value of L cep (5p) for 5p — > 0. The ratio R{Sp) then vanishes as 

R(Sp 0) = """I ( -^f f {tp)* (5.7) 



-2k B T cep A +t+ \pi- p v 

for 5p — > |^8| . The amplitudes of the power laws governing R(5p — > 0), which according 
to Eqs. ( |5T6| ) and ( |5.7| ) depend on the product /csT cep A a fe , show the same sensitivity to the 
type of the wetting agent (methanol-hexane or 4 He) as the effective Hamaker constant (see 
above). The drastic enhancement of kBT cep A a> b observed for typical binary liquid mixtures 
in comparison with 4 He makes critical effects on wetting layers much easier to detect exper- 
imentally. A corresponding statement can be made for direct force measurements by atomic 
force microscopes |j52l| . If two parallel plates at distance L are immersed into a binary liquid 
mixture, which is close to its critical demixing transition, the force per unit area K c between 
the plates will deviate from the bulk pressure due to the finite distance between the plates. 



This deviation is given by 48 



5K C (L) = K C {L) - K C (L = oo) = -— 5u(T c , L) = 2W e}f L^ (5.8) 



d_ 

if the van der Waals forces are not retarded [see Eqs. (|5.3| ) and (|5.4|) ]. Note that T c in Eq. 
( |5.8| ) is not given by T cep . Here T c marks a second order phase transition from the demixed 
to the mixed liquid, which takes place inside the liquid regime in the phase diagram away 
from the liquid- vapor coexistence surface (see Fig. 1 in Ref. j|9|])- However, typically T c is 
roughly about the same size as T cep . By inserting the values for A +i+ and A +i _ (see Table 
I), T c ~ 300i^, and W ~ 9 x 10 _15 erg for methanol-hexane into Eq. ( |5.4| ) one finds 

^ f -0.6 for (+, +) boundary cond. , , 

e ^ ' ~ 1 12 for (+, — ) boundary cond. 



According to Eq. (|5.9Q the critical contribution to 5K C (L) can lead to a sign reversal of 
5K C (L) for equal plates and increases SK C (L) by an order of magnitude for opposing plates. 
The effects of criticality on 5K C (L) should therefore be detectable by direct force measure- 
ments in critical binary liquid mixtures. 
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VI. SUMMARY AND DISCUSSION 



If macroscopic bodies are immersed in a critical fluid long-ranged forces between these 
bodies are generated by critical fluctuations of the order parameter. For the special case of 
binary liquid mixtures confined to a parallel plate geometry these forces have been analyzed 
for various boundary conditions involving surface fields in order to describe chemical affini- 
ties of the confining walls or interfaces towards one of the components of the mixture. In 
particular, the following results have been obtained: 

1. Within mean- field (Landau) theory for an Ising-like system (N — 1 order parameter 
components) the universal scaling functions f+,+ (y) and f+-(y) of the Casimir force can be 
easily obtained in a parameter representation without detailed knowledge about the order 
parameter profile. Either scaling function indicates that the corresponding Casimir forces 
should be visible over a surprisingly broad range in the scaling variable y = rl}l v . The 
scaling functions fsB,+(y) an d fo,+ (y) can be obtained from f +j+ (y) and f + _(y) by applying 
a simple scale transformation to / and y. In comparison with (+, +) and (+, — ) boundary 
conditions the Casimir forces for these mixed boundary conditions are substantially reduced 
both in their magnitude and in the range of the scaling argument y over which they are 
visible. For (+, +) and (SB, +) boundary conditions the force is attractive, for (+, — ) and 
(O, +) boundary conditions it is repulsive. For N > 2 an additional degree of freedom in 
the choice of the boundary conditions (surface fields) is provided by the introduction of an 
arbitrary tilt angle a between the surface fields. For N = 2 order parameter components 
and y = it is shown that the amplitude function g(a) smoothly interpolates between the 
special values (Casimir amplitudes) / +)+ (0) (a = 0) and / +i _(0) (a = n) of the scaling 
functions. The Casimir force vanishes for a = ir/3. For a = ir the order parameter profile is 
identical to the profile for N = 1 and (+, — ) boundary conditions. For critical binary liquid 
mixtures only the case N — 1 is relevant. 

2. For the special case y = (T = T Ct i, u ik) the scaling functions reduce to the universal 
Casimir amplitudes for (a, b) boundary conditions which have been calculated analyti- 
cally to one-loop order (Gaussian fluctuations) in order to obtain quantitative estimates for 
the magnitude of the Casimir force in d — 3. For the most relevant case N — 1 and for 
(+, +), (+, — ), and (O, +) boundary conditions it is possible to contruct Pade-type approx- 
imants for the Casimir amplitudes in d = 3 by including exact results from conformal field 
theory in d — 2 into an interpolation scheme for the amplitudes as a function of d. If a 3-loop 
estimate for the fixed point value u* of the renormalized coupling constant u is used in the 
interpolation scheme the resulting values for A + + , A +j _, and A 0j + in d — 3 agree quite well 
with corresponding numerical esitimates from a Monte-Carlo simulation of an Ising model 
confined to a slab geometry in d = 3 with surface fields. The estimates indicate, that for 
a critical binary liquid mixture the Casimir amplitudes are between one and two orders of 
magnitude larger than the previously studied amplitude Ao,o for 4 He at the A-transition. 

3. For critical binary liquid mixtures confined between equal or opposing walls the 
Casimir amplitudes A + + or A +j _, respectively, yield the absolute strength of the Casimir 
force in units of /cbT c . The film geometry considered here is realized in a natural way in the 
course of a wetting transition on a plane and chemically homogeneous substrate. The special 
case of (+, — ) boundary conditions is realized by the binary mixture methanol-hexane which 
forms a macroscopic wetting layer on Si - Si0 2 wafers in the vicinity of the critical end point 
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of the demixing transitions. Disregarding any temperature dependence of the Hamaker 
constant the presence of critical fluctuations in the wetting layer leads to an increase of the 
equilibrium layer thickness by more than a factor of two. The corresponding critical effect 
on a wetting layer of 4 He at the lower A-point is serveral orders of magnitude weaker. In 
accordance with this observation critical fluctuations in binary liquid mixtures have a strong 
effect on the effective Hamaker constant which determines the strength of the force between 
two parallel plates immersed into the mixture. Therefore, critical binary liquid mixtures 
appear to be ideal candidates to probe the universal Casimir amplitudes and the associated 
universal scaling functions by wetting experiments or by direct force measurements using a 
suitably adapted version of the atomic force microscope. 
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APPENDIX A: ORDER PARAMETER PROFILES 

The order parameter profiles in a critical film within mean field (Landau) theory for the 
Ginzburg - Landau Hamiltonian given by Eqs. (2.1) and ( j2.2|) have already been discussed 



in the literature in some detail for various reasons p. pT| . |32|j38[] (see also Sec. I). Therefore we 
only summarize the main results of mean field theory here for later reference. We restrict 
the analysis to the case N = 1 (Ising universality class). The Euler-Lagrange equation for 
the order parameter profile M(z) reads 

M"{z) = tM(z) + §M 3 (z), (Al) 

where the boundary conditions 

M'(0) = CiM(O) - h , M'(L) = -c 2 M(L) + h 2 (A2) 

must be fulfilled. In order to obtain the leading asymptotic behavior of M(z) in the critical 
regime we only consider the limiting cases h\ = h 2 —>■ oo [(+, +) boundary conditions] and 
hi = —h 2 — > oo [(+, — ) boundary conditions] in Eq. flX2]). In this limit the order parameter 
profile has the singularities M(z) ~ 1/z for z — > and M{z) ~ 1/ \L — z) for z — > L. This 
singularity of M(z) at the system boundaries just constitutes the mean field description of 
the asymptotic increase $(-s) ~ z~^l v of the order parameter profile as z — > for large (or 
infinite) surface fields. For this asymptotic power law to be valid the condition a<z<^ 
must be fulfilled, where a is a typical microscopic length scale and £ is the correlation length. 
In a lattice model for example a is given by the lattice constant. The order parameter profile 
for such a model will deviate from this power law increase on the scale z ~ a away from the 
surface and take a finite value right at the surface even for an infinite surface field. 

In order to simplify the notation for the following considerations we introduce the order 
parameter function m(z) by setting M(z) = \J\2/ g m(z) in Eq. ( |A1| ), where m{z) solves 
the modified Euler-Lagrange equation 
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m"(z) = rm(z) + 2m 3 (z). (A3) 
We furthermore suppress the parametric dependence of m(z) on the reduced temperature r 



in the notation. Multiplying Eq. (|A3|) by m'(z) one finds 



m' (z) = t m (z) + m (z) + m (zq) — r m (zq) — m (zq) (A4) 



as the first intgral of Eq. (|A3|) , where is an arbitrary reference point < z < L. 
For the combinations (+, +) and (+, — ) of boundary conditions considered here z Q = L/2 
is a convenient choice, because m(z) is either a symmetric or an antisymmetric function 
with respect to the midplane z = L/2, respectively (see also Refs. P^j2H|jr^] ). Up to an 
overall factor the integration constant in Eq. ( |A4| ) can be identified with (T±±) in the mean 
field approximation which we denote by (T±±) [see also Eq. (|3.2|)1. We define (Tj_j_)o = 
(6/(7) t±±, so that 

t±± = m' 2 (L/2) - t m 2 (L/2) - m 4 (L/2) (A5) 



is just the integration constant on the r.h.s. of Eq. (A4). With the substitution m 2 (z) = 
P(z) - r/3 Eq. (TAD takes the form 

P' 2 {z) = A[P{z) - e x ][P{z) - e 2 )[P(z) - e 3 ], (A6) 

where 

e x = -r/6 + y/T 2 /4-t ±± , e 2 = r/3 , e 3 = -r/6 - - t ±± . (A7) 



From the obvious property e x + e 2 + e% = and the structure of Eq. (|A6|) it is immediately 
clear that P(z) = m 2 (z) + r/3 is given by a Weierstrass elliptic function p(z; g 2 , gz) with 
the invariants 

g 2 = -4(eie 2 + e 2 e 3 + e 3 ei) = 4(r 2 /3 - t_Lj_), 

(73 = 4 ei e 2 e 3 = 4r(* ±1 - 2r79)/3. (A8) 

Moreover, p(z;g2,gs) has douple poles at z = and z = L, because m(z) has simple poles 
at these positions, so that the film thickness L is one of the periods of p(z; g 2 , #3). So far 
our statements are valid for both the (+, +) and the (+, — ) boundary condition. In order to 
derive the specific functional forms of the profiles we now consider each boundary condition 
separately. 

Turning to the (+, +) boundary condition first, we note that m'(L/2) = whereby 
t ±1 _ = -rm 2 (L/2) - m 4 (L/2) and Eq. flA7|) simplifies to 



e x = p{u t - g 2 , g 3 ) = P+,+(L/2) = m 2 +t+ (L/2) + r/3 

e 2 = p(uj 1 + u 2 ,;g 2 ,g3) = r/3 (A9) 
e 3 = P(u 2 ; g 2 , g 3 ) = -P+,+{L/2) - r/3 = -m 2 +>+ (L/2) - 2r/3. 

The quantities U\ and u 2 are the basic semiperiods of p[z\ g 2 , g 3 ). From Eqs. ( |A7| ) and 
(|A~9|) we conclude that e x > for all values of r and therefore P +>+ (z) = p(z;g 2 ,g 3 ) > for 
all < z < L. Therefore, the first basic semiperiod oj\ of the Weierstrass function can be 



23 



chosen as uo\ = L/2. It is then convenient to choose the second basic semiperiod U02 to be 
purely imaginary. We can now define the moduli k and k' of the corresponding Jacobian 
elliptic functions by |53| 



k 2 = = m )^\ + \ k- = l- k 2 . (A10) 

d - e 3 2m 2 (L/2) + r ' v ' 

According to Eq. ( |A10|) bulk criticality (r = 0) corresponds to k 2 = k' 2 = 1/2. The two 
basic semiperods are then given by the complete elliptic integrals of the first kind K = K(k) 
and K' = K(k') according to 



L K K .K> /A11 . 

ui = - = - = — , U) 2 = i-rrOJi. (All) 

2 y /e 1 - e 2 J2m 2 {L/2) + r K 



Combining Eqs. (|A10|) and (|A11|) we find the useful parameterization 



tL 2 = (2K) 2 (2k 2 - 1) , t ±± = -{2K/Lfk 2 {l - k ) (A12) 

of the Casimir force (Tj_j_)o = (Q/g)t±± as a function of the film thickness L and the scaling 
argument tL x I v = tL 2 within the mean field approximation. Finally, the order parameter 
function m +j+ (z) can be written in the form |53| 



2K dn(C; k) 2K 

m+ ^ z) = 17M^k) ' C = — z > (A13) 

where dn(£; k) and sn(^; k) are the Jacobian delta amplitude and sine amplitude functions, 
respectively. A slight disadvantage of Eqs. (|A12| ) and ( |A13|) is that in order to parameterize 
values tL 2 < —tt 2 one has to switch to negative values of k 2 , i.e., to purely imaginary moduli 
k in the Jacobian elliptic functions dn and sn. An alternative parameterization can be found 



easily by interchanging e 2 and e% in Eq. (|A7|) . From the corresponding modification of Eqs. 



( |A1(J| ) and ( |A11| ) we find the new parameterization (k 2 > 0) 

tL 2 = -(2K) 2 {k 2 + 1) , t ±± = (2K/Lfk 2 (A14) 

for rL 2 < — 7r 2 and the corresponding order parameter function reads 

2K 1 2K . . 

m +<+ W = -LMC^) ' C = ~L Z - (A15) 

From the symmetry of the order parameter profile for (+, +) boundary conditions it is 
obvious that within the mean field approximation the case of (SB, +) boundary conditions 
can be obtained from Eqs. ( |A12|) and ( |A13|) and their counterparts Eqs. ( |A14|) and ( |A15|) 



by the simple transformation L — > 2L. The corresponding order parameter profile is then 
given by m +i+ (z + L) evaluated in the interval < z < L. 

We now turn to the case of (+, — ) boundary conditions by noting that in this case 
m(L/2) = 0, because m(z) is antisymmetric around z = L/2. Therefore, we now have 



t±± = m (L/2) and instead of Eq. (|A9|) we find 
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d = p(ui; g 2 , g 3 ) = -r/6 - i\Jm+_(L/2) -t 2 /A 

e 2 = p(ui + u 2 , ; g 2 , g 3 ) = P+,-(L/2) = r / 3 

e 3 = p(<" 2 ; g 2l 9 3 ) = -r/6 + iyJm^_(L/2) - r 2 /4 



(A16) 



indicating that this time the two basic semiperiods are complex conjugates with u)\ + co> 2 
23R uji = L/2. In this case it is convenient to define the moduli k and k! as 



fc 2 = l/2-r/(4K +i _(L/2)|) , fc' 2 
The basic semiperiods can then be obtained from 



0J1 + uj 2 



L 

2 



L0 2 - L0 1 



1-k 2 



I — {UJl + co 2 ). 



(A17) 



(A18) 



y/\m' +t _(L/2)\ ' 

Combining Eqs. ( A17 ) and ( A18 ) as above we find the useful parameterization 

tL 2 = -2{2Kf{2k 2 - 1) , t ±± = (2K/Lf (A19) 

of the scaling argument tL 2 and the Casimir force (Tj_j_)o for (+, — ) boundary conditions. 
The corresponding order parameter function m+.-iz) can be written in the form [03 



2K cn(C; k) 
Xsn(C;A;)dn(C;fc) 



t 2K 



(A20) 



where in addition to Eq. ( |A13| ) the Jacobian cosine amplitude cn((; k) occurs. The param- 
eterizations given by Eqs. (|A19|) and (|A20|) have the disadvantage that values tL 2 > 2n 2 of 
the scaling variable correspond to purely imaginary values of the modulus k. However, in 
analogy with the (+, +) boundary conditions the alternative parameterization 



tL 2 = 2(2K) 2 (k 2 + 1) , t ±± = (2K/L)\l - k 



2\2 



(A21) 



can be found, where tL 2 > 2ir 2 corresponds to k 2 > and the corresponding expression for 
the profile m+-(z) reads 



2K cn(C;k)dn(C;k) 
L sn(C; k) 



2K 



-z. 



(A22) 



For (O, +) boundary conditions the Casimir force and the profile can be extracted from Eqs. 
( |A19|) and ( |A20|) or Eqs. (|A21|) and (|A22|) by the same simple transformation L — > 2L as 
described above for (SB, +) boundary conditions. 

We close this section with the remark that the order parameter profiles determined here 
can be written in the scaling form m(z) = L~^l v h{x\ y), where x = z/L and y = rL x l v are 
the scaling arguments and /3 — v — 1/2 within mean field theory. The y dependence of the 
profiles is determined by the above parameterizations y = y(k) in terms of the modulus k of 
the Jacobian elliptic functions. The scaling functions and can be easily read off 
from Eqs. (|A13|) and ( |A15|) and Eqs. (|A20|) and (|A22|) , respectively. 
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FIG. 6. Scaling function y) for T < T c (y < 0) (long dashed line), T = T c 

(y = 0) (solid line), and T > T c (y > 0) (short dashed line) according to Eq. ( A23| ) 
as a function of x. T c denotes the bulk critical temperature. For y ^ the thick-film 
limit (|y| 3> 1) is shown (see main text). 




FIG. 7. Scaling function h+-(x; y) for T < T c (y < 0) (long dashed line), T = T c 
(y = 0) (solid line), and T > T c (y > 0) (short dashed line) according to Eq. ( A24j ) 
as a function of x. T c denotes the bulk critical temperature. For y ^ the thick-film 
limit (\y\ 3> 1) is shown (see main text). 
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One obtains 



KA^y)= ^ n f Kx '^ , y=(2K) 2 (2k 2 -l), 
sn(2Kx; k) 



(A23) 



h +t+ (x;y)=2K \ , y = -(2K) 2 (k 2 + 1) 

sn(2Kx] k) 



and 



(A24) 



The functional forms of and below, at, and above bulk criticality are displayed 
in Figs. § and |7|, respectively. Bulk criticality means y = 0, i.e., A; 2 = 1/2 and off bulk 
criticality the thick film limit \y\ ^> 1 (k — > 1) is shown. In terms of the bulk correlation 
length £ the limit |y| 3> 1 in Figs.|] and [7] is represented as L/£ > 15. 

APPENDIX B: EIGENMODE SPECTRA 



The Gaussian Hamiltonian given by Eq. fl4.2|) can be conveniently diagonalized by solving 
the eigenvalue problem 

- V 2 ^(x, z) + 1(1 + l)m 2 (;z)^(x, 2) = £tf(x, (Bl) 

where / = 1 for the transverse spectrum and / = 2 for the longitudinal spectrum and 
< z < L. The film geometry is homogeneous and isotropic with respect to x so that we 
can write \I/(x, z) in the product form 

(x, z) = (27r)-^/ 2 e ip - x ^(z), (B2) 

where p is the longitudinal momentum and if)® (z) solves the eigenvalue equation 



+ 1(1 + l)m 2 (z)^(z) = egty?(*) (B3) 



so that the eigenvalue E in Eq. ( |B1|) takes the form E = p 2 + for I = 1 and / = 2, 
respectively. As shown in Eqs. ( |A6|) and (|A7|) m 2 (z) is given by the Weierstrass elliptic 
function p(z) = p(z;gz, #3), where g$ = for the case r = considered here [see Eq. 
( A~8D 1- Therefore Eq. ( P3|) is identical to the well known Lame differential equation f54l 



written in the form of an eigenvalue problem. The solutions of Eq. ( |B3| ) are known for 
I = 1 and I = 2 and can be used to construct the eigenf unctions ipjP(z). Note that due to 
m 2 (z) = p(z) ~ 1/z 2 for z — > one has ip$(z) ~ for z — > by inspection of Eq. ( |B3| ). 
Furthermore, (+, +) and (+, — ) boundary conditions can be treated on the same footing by 
noting that according to Eqs. ( |A11|) and ( |A18| ) one has 
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CO 



(+-) 



(1 - i)L/4 , 4 



(++) 



(+-) 



iL/2, 

(l + i)L/4 



(B4) 



for the basic semiperiods of the Weierstrass function. The spectra for the cases (SB, +) 
and (O, +) can be constructed from the spectra for (+, +) and (+, — ) boundary conditions, 
respectively. 

First we turn to the transverse spectrum . According to Ref. [54| the eigenf unctions 
ip^ (z) up to a normalization constant can be written in the form 



ij)V){z) = \a(z + a n )e~ zC(a " ) + a(z - a n )e z ^ an) ] /a(z) 



where 



-P(«n) 



(B5) 



(B6) 



yields the eigenvalues and ((z) and a(z) are the Weierstrass ( and a functions, respectively 
f53fl . The spectral parameter a n can be obtained from the requirement ip^(z) = ±ip^(z + 
L), i.e., the eigenfunctions are either even or odd functions when continued analytically to 
the interval [—L,L]. From Eq. |B|) one has L = 2uo{ ++) = 2(oj{ + ~ ] + c4 +-) ) and using; the 



shift properties of a(z) [|3|] the above shift operation can be directly applied to Eq. (|B5|) . 
One obtains for the eigenvalue spectrum 



2a n ((L/2) - L((a n ) = nici 



e (i) 



n > 2, 



where the lower bound on the mode index n comes from the requirement ip^ 
z — > for the transverse eigenfunctions (see above). 

al specti 



For the longitudinal spectrum (I = 2) the eigenfunctions take the form |54 



(2), 



where 



dz 



a n e 



+ a[z - a n )e 



z[Q(a n )+0„ 



p'[a r 



2p(a n ) + e ( ^/3 



p(a n ) 



=(2h3 



27<? 2 - 9(e 



(2)\2 



(B7) 
z 2 for 



(B8) 



(B9) 



yields the eigenvalues and p'(z) denotes the derivative of the Weierstrass p-function with 
respect to z. We again employ the symmetry requirement ipffl (z) = ±Vi (z + L) and the 
boundary behavior ipw (z) ~ z 3 for z — > to obtain 



2a n C(L/2) - L 



CK) + 



2p(a n )+el 2) /3 



nni , p(a n ) 



r(2)^3 



27^ 2 - 9(d z 



(2)x 



n > 3. (BIO) 



The solution of Eqs. ( [B~T| ) and ( |B10| ) for the eigenvalues eW, z = 1,2 cannot be obtained 
in a closed analytic form. In order to deal with the divergencies of the mode sums in Eqs. 
( |4.4| ) and ( f4.5| ) (see also Appendix C), we derive the asymptotic behavior of the eigenvalues 
from Eqs. ( |B7|) and ( [B10|) for large n. From the geometry of the problem it is clear that 
the leading term in an expansion of in powers of 1/n is given by the spectrum (rni/L) 2 
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of a free particle in a one-dimensional box of length L. Therefore the spectral parameter 
a n behaves as 1/n as w increases, so that the desired asymptotic form of the n dependence 
of the eigenvalues can be obtained from Eqs. ( P7| ) and ( |B10| ) by expanding the Weierstrass 
functions ((a n ), p(a n ), and p'(a n ) in powers of a n , where only the leading two terms are 
needed. Specifically, we use the expansions [^3| 

(Bll) 



C(x) = 1/x - g 2 /60 x 3 + 0(x 7 ) , p{x) = l/x 1 + g 2 /2Q x 1 + 0(x 6 ), 



where g% = is implicitly assumed. The calculation is straightforward so that we only briefly 
summarize the results for the eigenvalues eW. Corresponding expansions are obtained for 
the spectral parameter a n which will not be reproduced here. 
For (+, +) boundary conditions one has 



C(L/2) = CK 



(++h 



tt/(2L) , g 2 = -At ±± = (2K/LY 



(B12) 



where K = K (1/^/2) [see Eqs. (|X8|), and Q)]. By insertion of Eqs. (|BTl| ) and ( BT2l) 

into Eqs. ( P37|) and (|B10|) one obtains the expansions 
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For 



boundary conditions one has correspondingly 



aL/2) = a^ ) + 4 + - ) )=n/L, g 2 



-At 



-A(2K/LY 



where K is given as above [see Eqs. (|7^8|), ( |3.3| ), and (|3.5|) 1. Insertion of Eqs. 
( P14| ) into Eqs. ( p7|) and ( |B10|) yields the expansions 
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The asymptotic expressions for the spectrum given by Eqs. ( B13|) and (|B15| ) capture all 
divergent terms in the mode sums in Eqs. (4^) and ( [4.5|) as will be seen in Appendix C. 
Furthermore, Eqs. (pi3| ) and ( pi5|) provide very good initial values for a numerical solution 



of Eqs. (|B"T| ) and ( |B10| ) by iterative schemes, e.g., the Newton procedure. 

For (SB, +) boundary conditions the eigenvalue spectra can be obtained from the case 
of (+, +) boundary conditions by employing the transformation L — > 2L and by allowing 
only even indices n for e^ 1 ) and only odd indices n for [see Eq. (|C7|) 1. Likewise, the 
eigenvalue spectra for (O, +) boundary conditions can be obtained from the case of (+, — ) 
boundary conditions by again employing the transformation L — > 2L and by allowing only 
odd indices n for and only even indices n for [see Eq. flClO )]. The reason for this 
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simple rule is that for (+, +) boundary conditions starting from the ground state every 
second eigenfunction has vanishing slope at z = L/2 so that after rescaling L — > 2L the 
eigenfunctions for (SB, +) boundary conditions are already contained in the (+, +) case. An 
analogous argument relates the spectra for (+, — ) and (O, +) boundary conditions starting 
from the first excited state for the (+, — ) case. 



APPENDIX C: REGULARIZED MODE SUMS 



The mode sums appearing in Eqs. ( |4.4j ) and ( [4.5[ ) are divergent for any spatial dimension 
d of interest. Within the dimensional regularization scheme used throughout this investi- 
gation d is used as a free parameter in order to find an analytic continuation of the mode 
sums as a function of d, where d = 4 — e is this case. On the other hand the mode sums in 
Eqs. ( |4.4| ) and ( [4.5| ) also constitute the zeta functions of the eigenvalue spectrum with a d 
dependent argument [55j. The zeta function regularization of mode sums, which is a widely 
used technique to treat divergent series like those in Eqs. (|4.4| ) and Q4.5Q |55| , is therefore 
equivalent to the dimensional regularization scheme. 

The major obstacle towards an analytical treatment of the aforementioned mode sums 
has been removed in Appendix B by the derivation of the asymptotic behavior of the eigen- 
value spectrum for large mode numbers given by Eqs. ( B13| ) and ( P15| ). Using these results 
one has for i — 1,2 




(ih(d-i)/2 



'nn' 


d-l 


1 - 


2A 


B- 




[ L _ 






h n 1 . 


(d-l)/2 j 



E 

n=no 



n 



d-7 



(CI) 



which is convergent for any d of physical interest and can thus be determined numerically 
from the solutions of Eqs. ( |B7| ) and ( |B1(J| ) for the transverse and the longitudinal mode 
sum, respectively. The problem of regularizing the mode sums has therefore reduced to the 
regularization of the corresponding sums over the large n expansions given by Eqs. ( |B13| ) 
and ( |B15| ), i.e., one has to consider the series 
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for d — 4 — e. If the lower summation bound n in Eq. flC2j ) is chosen sufficiently large, one 
can safely expand the term under the sum in powers of l/n 2 which leads to an expansion of 
the series given by Eq. ( |C2| ) in terms of Hurwitz functions ((x, no). One finds for d = 4 — e 
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where the e-expansion has already been carried out up to terms 0(e). The expansion shown 
in Eq. ( |C3|) converges quite fast already for 3 < n < 5. The 1/e pole indicating the UV 
singularity can be extracted from Eq. (|C3|) using the expansion 



no-l 



C(l + e, n ) = 1/e + 7 - E 1 A + O(e) 

fe=i 



(C4) 



where 7 ~ 0.577216 is the Euler constant and Uq is a positive integer. With the coefficients 
A and B taken from Eqs. ( P13|) and ( |B15|) the expressions given by Eqs. ( Ol]) , 



and 

can be combined to the following regularized and e-expanded expressions for the mode 



sums. 



For (+, +) boundary conditions one finds with K = K(l/y/2) 
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0(e). 



For (+, — ) boundary conditions the corresponding result reads 
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For (55, +) boundary conditions we apply the simple transformation described in the last 
paragraph of Appendix B to the eigenvalue spectrum for (+, +) boundary conditions. From 
Eq. ( |B13| ) we find the expansions 
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(n>l), 



for the transverse and the longitudinal spectrum, respectively. Due to the appearance of 
half-integer arguments in the Hurwitz functions for the transverse mode sum in this case 
one needs the expansion 



C(l + e,3/2) = l/£ + 7 + 2hi2-2 + C>(£) 



(C8) 



instead of Eq. (|C4|) . Furthermore, the r.h.s. of Eq. (|C3|) with n replaced by 3/2 is 
needed in order to derive the regularized longitudinal mode sum. The transverse mode 
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sum, however, can be evaluated directly using Eqs. ( |C3j ) and ( |C4|) . One therefore finds for 
(SB, +) boundary conditions 
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For (O, +) boundary conditions we apply the same transformation to the eigenvalue spec- 
trum for (+, — ) boundary conditions. From Eq. (P15|) we find the expansions 
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Using Eq. ( |C3|) with no replaced by 3/2 in order to evaluate the transverse mode sum one 
therefore finds for (O, +) boundary conditions 
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In order to facilitate the e-expansion of Eq. ( |4.5| ) we finally note that 
T [(3 - d)/2] = -2^ [l + £ (1 - In 2 - 7 /2) + C(e 2 

and 

r(d/2) = i - £-(i - 7 )/2 + o(c- 2 ) 

for d = 4 — e (see also Ref. |28j for similar relations). 
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